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We investigate the dynamics of a one-dimensional asymmetric exclusion process with Langmuir 
kinetics and a fluctuating wall. At the left boundary, particles are injected onto the lattice; from 
there, the particles hop to the right. Along the lattice, particles can adsorb or desorb, and the 
right boundary is defined by a wall particle. The confining wall particle has intrinsic forward and 
backward hopping, a net leftward drift, and cannot desorb. Performing Monte Carlo simulations and 
using a moving-frame finite segment approach coupled to mean field theory, we find the parameter 
regimes in which the wall acquires a steady state position. In other regimes, the wall will either 
drift to the left and fall off the lattice at the injection site, or drift indefinitely to the right. Our 
results are discussed in the context of non-equilibrium phases of the system, fluctuating boundary 
layers, and particle densities in the lab frame versus the frame of the fluctuating wall. 

PACS numbers: 



I. INTRODUCTION 

Asymmetric exclusion models with a fixed [H,[l,[l] , 
and typically large number of lattice sites have been 
the subject of much recent theoretical attention 
|, i, I 0, 1 1 E3, El E3, H 0, El- Biophysi- 
cal applications and new fundamental understanding 
of non-equilibrium steady-states (NESS) have moti- 
vated many extensions of the simple totally asym- 
metric exclusion process (TASEP) with open bound- 
aries. These include partially asymmetric models, 
where particles can hop backward [H , exclusion pro- 
cesses with nonuniform hopping rates [3, [3, 0"3] , ex- 
clusion among particles of arbitrary size Tfl, 3, H3] > 
multispecies exclusion processes [3, E3, Gj],ll2|, mul- 
tichannel exclusion processes [l8j] and exclusion pro- 
cesses with L ang muir t y pe adsorption and desorp- 
tion kinetics [l3l [Til Il5j . All of these studies have 
considered open, well-defined boundaries, where the 
length of the lattice is fixed. TASEP models with 
one open and one closed boundary conditions have 
also been considered [Hj] ■ 

However, applications may arise where the length 
of the system is dynamically varying. The sys- 
tem size may vary because a single particle pushes 
against a boundary-defining wall. One example is 
helicase- induced opening of replication forks in DNA 
processing 20]. Here, the moving replication fork 
defines a moving boundary of the system. Examples 
of variable-system size exclusion processes that in- 
volve multiple motor particles include mRNA trans- 
lation in the presence of hairpins in the mRNA, and 
molecular motors processing on elongating actin fil- 
aments. Ribosomes that process along mRNA (in 
the process of protein synthesis) during translation 
[3, HH often encounter a hairpin and the position 
at which the hairpin starts represents a wall over 



which the processing ribosomes cannot pass. The 
detachment rates of the ribosomes and the tightness 
of the hairpin may determine if the ribosomes can 
translate the mRNA through the hairpin sequences. 
Actin polymerization at the leading edge of filopodia 
also seems to be mediated by processing molecular 
motors that may carry actin assembly components 
[22l . [23l | . The motors detach, and possibl y a ttach, 
anywhere along the growing actin filament [24 ] . The 
depolymerization of the leading tip may be limited 
or enhanced by the presence of a motor or other 
actin associated proteins [25l l2t| . Finally, a recent 
model of a dynamically extending exclusion process 
without Langmiur kinetics has been recently stud- 
ied [13] • This model has been applied to filamentous 
hyphae growth in fungi [28j . 

With the above applications in mind, we consider 
a TASEP with a dynamically varying length. Specif- 
ically, we analyze a many-particle asymmetric exclu- 
sion process with a fixed open boundary on the left, 
a fluctuating boundary on the right, and Langmuir 
kinetics. The particles have a fixed injection site 
and can adsorb and desorb. A wall with an intrinsic 
leftward drift (representing e.g., a hairpin which en- 
ergetically favors spontaneous closing or the barbed 
end of an actin filament that prefers depolymeriza- 
tion) prevents the passage of particles. The particles 
advance and provide a pressure against the wall. For 
certain attachment/detachment and wall hopping 
rates, the system reaches a NESS in which the statis- 
tics of the wall position are stationary. For other val- 
ues of the kinetic parameters, no time-independent 
mean wall position exists. The wall will either drift 
steadily towards the particle injection site and fall 
off the lattice, or move indefinitely away from the 
injection site, continuously increasing the size of the 
system. The specific details of the stochastic pro- 
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cess are shown in Figure Q] Particles are injected 
into the first lattice site with rate a provided it is 
empty. In the interior of the lattice, each particle 
moves forward with rate p only if the site ahead of it 
is unoccupied. Particle attachment and detachment 
occur with rate k + and fc_ , respectively, throughout 
the lattice. 

a k_ k. 
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FIG. 1: A totally asymmetric exclusion process bounded 
by a fluctuating wall. Particles are injected onto the 
leftmost site with rate a, and move to the right with rate 
p. In the interior, particles detach and adsorb with rates 
fc_ and respectively where k± <C p- The lattice is 
bounded on the right by a fluctuating wall with intrinsic 
hopping rates w+ and w-, where w+ < w-. 

The lattice length is not fixed, and TV denotes the 
position of the particle-confining wall that hops for- 
ward with rate w+ , and backward with rate W- pro- 
vided there is no particle to its immediate left. The 
particle occupation at each site, 1 < i < N — 1, left 
of the wall is represented by the occupation vari- 
able Oi £ {0, 1}. If w_ < w+, the wall will move 
indefinitely away from the injection site. In order 
to prevent the wall from always escaping to infinity, 
we consider the more interesting case of an intrinsic 
leftward drift described by iu_ > w + . 

The wall position TV is not fixed (even at steady 
state), but rather, is determined by the intrinsic wall 
hopping rates, and the exclusionary interactions be- 
tween the wall and the lattice particles. Our analy- 
sis is aimed at understanding how the wall dynam- 
ics depend on the parameters a,k±,w±,p. In the 
next section, we derive relations for the distribution 
functions of the wall. In steady-state, these relations 
constrain the particle density at the wall. In Section 
IIII1 we use mean field theory (MFT) to solve for 
the density profile, and show that the density profile 
obtained using mean field theory is inaccurate near 
the wall. Since quantitative prediction of the wall 
dynamics will require accurate determination of the 
particle densities near the wall, in Section IIV1 we 
develop a moving-frame finite segment mean field 
approach to accurately solve for the density profile 
near the wall. The existence of a steady-state solu- 
tion and the dependencies of the mean wall position 
(TV) on the problem parameters are explored and 
plotted in Section IVl 

II. WALL DYNAMICS 

The net drift of the wall is the difference be- 
tween its forward and effective backward hopping 



rates. The effective backward hopping rate depends 
on both the intrinsic backward hopping rate u>_ , and 
on the occupancy of the site immediately to the left 
of the wall since a particle there will block the wall's 
backward motion. The wall's rightward hopping is 
never impeded. The probability of finding a particle 
immediately to the wall's left varies with its position, 
thus, the wall dynamics are position-dependent. De- 
fine QN(t) as the probability that the wall is at posi- 
tion TV at time t, and Q' N {t) as the joint probability 
that the wall is at position TV at time t and the site 
just before the wall is empty. The wall dynamics 
obey 



Q N {t) = w-Q' N+1 - w^Q' N - w + Q N + w + Q N ^x, 

and the moments of the wall position can be formally 
expressed as 

i 00 fc— l /u\ 

-(N k ) = -£(w^Q' N -w + Q N )J2(-V k - j ( ;)N j . 

(2) 

Although one cannot find Qn or Q' N explicitly with- 
out solving the full exclusion problem, we can take 
k = 1 in @ to determine the mean wall velocity via 

-(N)=-W-J2Q'n + w+. (3) 

' JV=0 

If the mean wall position is time-independent, 
S5v=o Q'n = w +l w -i an d the expected occupancy 
of the site immediately preceding the wall is 

oo 

( (TJV _ 1 ) = 1-^Q^ = 1-H±. (4) 

We show in section IV Bl that there are some pa- 
rameter regimes in which Q cannot be satisfied. 
For these parameter values, there exists no time- 
independent mean wall position. However, one can 
still use ([31) to determine the relevant mean wall dy- 
namics. The preceding analysis suggests that it may 
be more natural to define sites near the wall by their 
position relative to the wall than by their absolute 
position on the lattice. To avoid working in both 
frames of reference, in the next section, we will begin 
by considering the limit in which the wall hopping 
rates are small compared to other rates in the prob- 
lem (p, k + , and fc_). In this limit, the wall dynamics 
are slow compared to the particle dynamics, and we 
will assume that the wall frame is stationary. 
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III. MEAN FIELD SOLUTION OF DENSITY 
PROFILE 

In the w±/k±, w±/p — > limit, we expect the wall 
to be nearly stationary. Mean field equations can 
be derived by ensemble-averaging the rate equations 
for the occupation variables <jj, and ignoring corre- 
lations ((<7j<7j) ~ (<7i)((jj)). Upon defining the mean 
occupation Si = (<7i), the mean field equations for a 
fixed (w± = 0) wall system in NESS are 



— j = -Si(l - Si+i) + Sj_i(l - Si) - k-Si 
+fc+(l - = 0, 



ds 

— = a(l - sx) - fc_s! - si(l - s 2 ) 
di 



(5) 



di 



+fe f (l-si) = 0, (6) 

-fc_SjV-l + - Sjy-l) 

+«w-a(l - sjv-i) = 0. (7) 



where the adsorption, desorption and injection rates 
have been normalized by p and time has been 
rescaled by p _1 - hence, k±, a and t in (07]) are 
dimcnsionless. 

However, in order to use condition ([4]), we need 
expressions for particle density at sites defined by 
their distance from the wall. In the fluctuating frame 



of the wall, we use the notation s, 



SN-i 



Upon 



rewriting ([5][7]) in the wall frame, we find 



dt 



ds N . 



dt 

dsi 



-(1 + W + )Sj(l - Sj-i) + (1 + (1 - Si)uu)Sj+i(l - Sj) - k-8j + fc+(l - Sj) 

-w-(l - Si)(l - s j+ i)sj + w+(l - Sj)sj-i = 0, (8) 

a(l - s N -i) - /c-sat-i - sjv-i(1 - sjv-2) + fe+(l - sat_i) = 0, (9) 

-k-h + k+(l - h) + 0- + w-)h(l - h) - hw+ = 0. (10) 



As expected, {8j and {TO} reduce to © and in 
the w± = limit. If the position of the wall were 
fixed, we could simply use the iteration given by ([5]) , 
along with boundary conditions © and {7} to solve 
for the density profile s^. 

Now consider a moving wall problem. Because 
(N) is undetermined, we need three conditions to 
solve ((HJ. In addition to the two boundary con- 
ditions ([9]) and (|10|) . we require a third condition, 
Si = l — w+/w—, to determine (N). This third 
boundary condition fixes si; §2 is set by (|10l) . and 
we can use ([8]) to iterate forward in j as many times 
as required toward the injection site, until ([9]) is sat- 
isfied. The number of iterations required to satisfy 
© determines the mean position, (N), of the left 
boundary, and hence the NESS size reached by the 
system. Although ([5]) was derived in the wall frame, 
the resulting density profile is nearly identical to a 
stationary frame profile derived from ([5|) when Sj is 
not varying rapidly with site i. See the Appendix 
for further discussion. 

For standard particle-conserving TASEP models, 
away from boundaries, MFT predicts the particle 



densities to a very high accuracy [2| • A conservation 
law for the particle density can be used to fix the 
end densities to their exact values so that the MFT 
also performs well near boundaries [lf|. In Figure 
UK, we plot the density profiles from Monte-Carlo 
simulations and mean field recursion relations for 
the simple TASEP (A;± = 0) with a fixed number of 
sites, N — 10000. Differences in the density profiles 
are evident in the insets. 

Because we include particle adsorption and des- 
orption through Langmuir kinetics, there is no con- 
servation law for the particle density. In this case, 
the boundary densities are not fixed and we see in 
Fig. that simple mean field calculations of the 
boundary density can differ appreciably from the 
values found from Monte-Carlo simulations. How- 
ever, MFT still matches simulation results in the 
bulk where s$ varies slowly. In the following section, 
we use an approach that couples explicit enumera- 
tion within a finite segment of sites to the mean field 
results accurate outside the segment. This finite seg- 
ment mean field theory (FSMFT) includes particle 
correlations within a segment of sites adjacent to the 
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FIG. 2: (Color online) A comparison of density pro- 
files derived from Monte-Carlo simulation and MFT. (a) 
The MFT and MC density profiles for conserved particle 
TASEP are compared (a = 0.6, (3 = 0.9, N = 10,000). 
Despite differences between MFT and MC in the bound- 
ary layers, the particle density at the ends (i — 1, 
i = 10000) are matched through particle conservation. 
Insets show the left and right boundary layers in de- 
tail, (b) For a TASEP with Langmuir kinetics (with 
k- = 0.01, N = 300, p = 1, and k+,w+,W- = 0), the 
MFT density profile can be appreciably different from 
the MC results, especially near the boundaries. 



IV. FINITE SEGMENT METHOD 



We have shown that mean field theory does a 
poor job of predicting the profile Si~jv-i near the 
wall when there is a boundary layer. To more accu- 
rately compute the particle density in this region, we 
will solve the Master equation for a finite segment 
of m sites preceding the wall. First, we introduce 
some notation to explain the mechanics of the finite- 
segment mean field theory (FSMFT). For the binary 
string (<7jv_ m , cn-i), corresponding to the 

occupancy of sites in the finite segment we define 
the state of the segment as the base ten value of the 
string. For example, for m = 2 sites just left of the 
wall, we have four possible combinations for the oc- 
cupancies (00), (01), (10), and (11) corresponding to 
states i = 0, 1, 2, 3, respectively. If Pi is the proba- 
bility that the finite segment configuration is in state 
i, the Master equation is dtPi — MijPj where My 
is the transition matrix. In the m = 2 case, 



M 



/ -(1 + w-)s* - 2k+ k- 

k + — (w + + fc_ + k^ 
(1 + w_)s* + k + w + 
\ k+ + s* 



k- +w+ \ 

s* l + w-(l-s*) k_ 

— {w + + W- + 1 + fc_ + k + ) k- + w + 

k+ + w-s* —2k- — u>+ / 



where s* = (<JN-m-i) is the mean occupancy in the 
lattice site just to the left of the segment. The mean 
occupancies in the finite segment can be calculated 
from M in the following way. First, the eigenvec- 
tor, p(°), corresponding to the eigenvalue zero is 
computed. The vector P", normalized such that 



(11) 



r 



S™ 1 = 1 corresponds to the stationary prob- 
ability distribution, i.e., (9 t P(°) = 0. Let v be a 
m x 2 m matrix where the columns are the ordered 
state vectors. The mean densities are then given by 
(sjv-to-i, • ■ • , sjv-2, s W -i) T = vP (0) . 



5 




0.01 0.02 0.03 0.04 0.05 0.06 
S* 



(b) 




- - Mean Field Theory 


I 

1 


— FSMFT 


1 


♦ Monte Carlo 


1 

1 ' 
1 / 
1/ 
1/ 

f\ 
\ 


A 


4 1 

/ 1 
/ 1 
/ 1 
/ 1 



292 294 296 298 300 

position i 

FIG. 3: (a) F(s*) is plotted with parameter values 
k+ = 0, k- = .01, w+ = 0.005, W— = 0.01. (b) The 
finite segment method predicts the boundary layer pro- 
file significantly better than mean field theory. Here, the 
final ten sites of a fixed-wall profile (TV = 300) are plot- 
ted. The parameters k- = 0.01, k+ = 0, and a = 1 were 
used. 

For every value of s* , FSMFT can be used to 
compute the mean densities sjv- m , sjv-2> s/v-i- 
In particular, it establishes a one-to-one correspon- 
dence between s* and Sjv-i: 

SJV _! = F(s*;w±,k±). (12) 

Our calculations indicate that F(s*) is always a 
monotonically increasing function of s*, as shown 
in Figure Comparing the density profiles near 
a fixed wall (u>+ = u>_ = 0), Figure [3Jd shows that 
using FSMFT (with m = 5) significantly improves 
our prediction of the particle density near the wall 
over that obtained using simple (m = 1) MFT. To 
calculate (N) for a fluctuating wall (w- > w+ > 0), 
we first solve for the profile of a segment of sites ad- 
jacent to the wall. From when the wall attains 
a steady-state position, sjv_i = 1 — w + /w~. Using 
(TT2")) , we find the value of s* satisfying 1 — w+ / w_ = 
F(s*). Defining this particular value of s* as s* 9 , we 
then use the recursion relation given by {8| to solve 
the density profile to the left of the finite segment. 
Using the values of SN- m and s* eq = SN- m -i from 



the FSMFT as starting conditions for the recursion 
equations, we iterate to the left until the left bound- 
ary condition ^ is satisfied. In summary, the finite 
segment mean field theory (FSMFT) is implemented 
by the following steps: 

• For a given s*, solve for the normalized eigen- 
vector corresponding to the zero eigenvalue of the 
2 m x 2 m transition matrix M tJ (s*). 

• From the zero eigenvector, express the mean 
density sn-i at the site nearest the wall as a func- 
tion of s*, giving relation (TT2")) . 

• For a static wall NESS, set sn-i = 1 — w + /w- 
and find s* q that yields zero net wall drift by using 
1 - w + /w_ = F(s* eq ; w±,k±). 

• Starting with s m+ i = s* eq (and s m = s N _ m ) 
iterate using the simple mean field equation ([S]) until 
equation © is satisfied. 

• The number of iterations required determines 
the mean wall position ((N) ~ number of iterations 
+to+2) as a function of the rate parameters through 
the starting value s* eq . 

We expect the predicted results from a moving- 
frame FSMFT to be in good agreement with those 
from MC simulations. This is because in regions 
where s is slowly varying, the mean field equations 
describing the density profile in the wall frame and 
in the lab frame yield nearly identical profiles. This 
can be seen from the continuum equations, as will 
be discussed in the Appendix. In these regions, ac- 
curate estimates of the mean wall position can also 
be obtained using the continuum approximations to 
©, provided (N) is large. When state enumera- 
tion of a larger segment is used, more of the cor- 
relations within the density boundary layer is taken 
into account and more accurate results are expected. 
Provided most of the regions with large gradients in 
density is captured by the finite segment, the re- 
sults will be very accurate. The incremental accu- 
racy achieved as larger segments are used has been 
discussed in a different, but related system @. In 
the subsequent analyses, we use a five-site (m = 5) 
FSMFT - generating a 2 5 x 2 5 eigenvalue problem 
in the process - and self-consistently solve for the 
densities away from the boundary layer. This choice 
of segment size is sufficient to yield accurate results 
for all parameters explored. 
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V. RESULTS AND DISCUSSION 

A. Time-Independent Mean Wall Positions 

We first consider regimes in which the wall ac- 
quires a static mean position in NESS. Using Monte 
Carlo simulations and FSMFT, we study the depen- 
dence of the mean wall position on the injection rate 
a, particle adsorption and desorption rates k±, and 
the wall hopping rates w±. We can use analytic 
solutions of the bulk continuum equations in order 
to understand parameter dependencies of our model. 
Although mean field theory poorly describes our sys- 
tem in the boundary layers where the profile varies 
rapidly, away from boundary layers, simple MFT is 
accurate. In these regions, to guide our our anal- 
ysis, we will use the continuum limit of the mean 
field equations. We define e = 1/Nq where No is a 
characteristic number of lattice sites (to be derived 
below) and x = (i — 1) /No as a relative position 
along the lattice. As shown in the Appendix, the 
NESS density profile obeys 
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e{2s - l)s'(x) +k+- (k+ + k-)s + 0{e 2 ) = (13) 

in both the lab and wall frames of reference. Upon 
integrating, we obtain the implicit equation 

(fc+ - fe_) In \k+ - (k+ + k-)s\ - 2fc+ + 2(fc+ + k-)s 



{k+ + k-f 



C, 



(14) 

where C is a constant of integration. In the con- 
tinuum description, the entrance site is at position 
x = 0, the wall's position is L, and the mean wall 
position is (L) = e(N). We can use (Tl4"|) to under- 
stand the behavior of the mean wall position (N). 
First, note that the left hand side of fT4]) scales as 
(k + + /c_) _1 . Since e _1 = N scales as (fc_ + fc+)~\ 
we define Nq = {k + + fc_) _1 . For a continuum de- 
scription to be useful, Nq must be large, so k± must 
be small. 

Equation (|14p gives an implicit formula for the 
bulk density, which we denote by Sb{x), in terms 
of the adsorption and desorption rates, k± and the 
integration constant C . The injection rate a deter- 
mines C, and along with the wall hopping rates w±, 
determines the mean wall position. As shown in the 
Appendix, the solution near the left boundary varies 
slowly when a < 0.5. Furthermore, if fc_ <C a, we 
can approximate s\ « S2 in J5[ to conclude that 
s(0) w a. This simplified condition can be used to 
determine C in (| 14|) . 

When a > 0.5, a boundary layer arises on the left 
(this can be seen in Fig. Sh.). In this regime, s(0) 



FIG. 4: (Color online) Simulations were performed with 
w+ = 0.001, W- = 0.01, k- = 0.01, and different values 
of a. In (a), profiles for three simulations are plotted 
with their mean wall positions aligned. When a < 0.5, 
s(0) w a. When a > 0.5, ss becomes multivalued and 
there is a boundary layer on the left. Within the bound- 
ary layer, a small change in the position results in a large 
change in the particle density. Thus, in (b), where (iV) 
is plotted as a function of a, large changes in a result 
in small changes in (N) when a > 0.5. Also shown is 
the prediction from (|16|l with s* eq = 0.026, determined 
by FSMFT. 



can no longer be approximated as a, and sb(x) be- 
comes invalid near the injection site. While sb(x) 
is still a good approximation to the density profile 
outside the boundary layer (where s ^ 0.5), there 
is no straightforward, analytic way to calculate C 
when a > 0.5. In Fig. [4ji, when a = 1, C is used 
as a single fitting parameter and is determined em- 
pirically such that in the bulk, sb(x) approximates 
the density profile obtained using MC simulations. 
The mean wall position, (L), is found through the 
relation sb{(L) — (m + 1)) = s* eq where s* eq is found 
using an m— site FSMFT and is the value of s* that 
puts no net drift on the wall. 

Figure 2k shows results from MC simulations, 
shifted so that the mean wall positions are aligned 
at (N) — 225, which is the mean wall position when 
a = 1. While the density profile has a sharp bound- 
ary layer at the wall in the wall frame, in the lab 
frame, the boundary layer is smeared out due to wall 
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fluctuations. This results in the broad peaks cen- 
tered on the mean wall position seen in Figure [4^,. 
The outer solution Ss(i/No) with N = (fc+ + fc_)~ 1 
is shown by the dotted curve. The close agreement 
between the MC data and sb (x) suggests that drop- 
ping the 0(e 2 ) term in (fT3|) to obtain sb(x) produces 
an excellent approximation to the mean particle den- 
sity, provided a < 0.5. Note that a, through C, sim- 
ply shifts sb(x) to the left or right; thus, when we 
vary only a and plot the resulting density with the 
mean wall positions aligned (as they are in0k), the 
profiles collapse onto the same curve. 

We can also use (TH|) to predict the mean wall 
position as a function of the injection rate a when 
a is not too large. For simplicity, consider k + = 
- the analysis for k + ^ is analogous. Using the 
simplified condition s(0) — a in (HU, we have 



C = 



2a — ln(afc_ 



Now, using the relation s* eq = sb{(L) — e(m 
and ([TBI . (fT4| becomes 



(A0 = r-ln 



a e 



2s!„ 



(15) 
M)) 



(16) 



The dependence of (N) on a is shown in Figure 
|4)d, predicted using four different methods. Simple 
MFT (m = 1, dotted curve) performs poorly rela- 
tive to MC simulations (open diamonds). The re- 
sults from FSMFT with m — 5 (solid curve) agree 
very well with the MC data for all values of a. The 
solution of (fl6|) (dashed curve) performs reasonably 
well provided a is not too large. When a > 0.5, 
s(0) = a is a poor approximation to ([6]) and the re- 
sulting prediction of (N) suffers. In fact, the slope 
s' B diverges when sb = 0.5, which can be seen from 
(flU)) . When a > 0.5, there is a boundary layer on the 
left with width O (y/e) (cf. Appendix). As a result, 
increases in a above 0.5 will increase the height of 
the boundary layer, but will not significantly change 
the mean wall position, and (N) becomes insensitive 
to changes in a [3l| . 

We now discuss how changes in the wall hopping 
rates can affect the wall position. In Fig. [SJi, for 
a fixed value of w+, one sees that an increase in 
w_ increases the value of s* eq . Our FSMFT pre- 
dicts that given values of w+ and , s* eq must sat- 
isfy 1 — w+/w- = F(s* eq ]w±). For small values of 
w±, F(s+ q ,w±) w F(s* ? ,0) Ml- suggest- 
ing that s* q depends primarily on the ratio w + / w_ , 
with only a weak dependence on the individual wall 
hopping rates. Since F is a monotonically increasing 
function, s* eq increases with u>- / w+ . 

A change in s* eq induces a change in the mean 
wall position, shown in Fig. [5Jd. Again, this is 



consistent with our theory since (L) must satisfy 
sb((L)— e(m+l)) = s* g . In the special case a < 0.5, 
one can use (fl6| to predict (N) directly, given s* g . 
When a > 0.5, one either has to solve the full set of 
discrete MFT equations ((8j) and (O - or the equiv- 
alent continuum equations l|28p and (f!?5)) - coupled 
to a finite segment, to obtain (N). Our results from 
solving the discrete equations are shown in Fig. [5l[c) . 
Using simple MFT (m — 1) without a larger fi- 
nite segment generally results in poor predictions for 
(N). 

A more complete understanding of the wall dy- 
namics can be garnered by analyzing the wall fluc- 
tuations. For simplicity, we consider the contin- 
uum description in which the wall's motion can 
be approximately described by a diffusion constant 
D = e 2 w + and a position-dependent drift, V(L). If 
one assumes that the wall fluctuates within a har- 
monic "potential," this drift takes the form V(L) — 
-a(L - (L)), where a = -(dV/dL)\ {L) . This ap- 
proximation effectively closes (TTJ) by expressing the 
effects of conditional probability, Q' N , in terms of 
a drift. For L w (L), the probability density of the 
wall's position, Q(L) (the continuum analog of Qn), 
can be approximately found from the solution of 



dQ{L,t) 
dt 



(L))Q(L)]+D 



d 2 Q 
dL 2 ' 



(17) 



Upon imposing the normalization J_ Q(L)dL = 1, 
we find the steady-state solution to (jTTJ) : 



Q(L) = 



where a is given by 



e 2D 



2Dn 



(18) 



dV_ 

111 



BV f)F 

= -Wd? s ' B{{L) - £{m+l)) - (19) 



The drift velocity V(L) can be inferred from 
e(w + — iu_(l — sjy-i)) and p^|) . which relates the 
mean occupancies at positions L and L — s(m +1), 



V(L) 



e idj 



v- + w-F(s*(L - e(m + 1)))] . 

(20) 

By defining the drift V(L) using the steady-state 
relation F, we have implicitly made an adiabatic 
approximation where the particles have reached a 
NESS for any wall position L. 

We can also estimate the variance of the wall po- 
sition by using E 2 w D/a, and (TT5)) for a. Upon dif- 
ferentiating V(L) — e(w + — W- [1 — F{s' B {L — e(m + 
1)]), we find (dV/dF) = ew-. We can estimate 
(dF/ds*)\ s * using the finite segment method, and 
we know s' B ((L) — e(m + 1)) exactly from (fT3|) . As- 
suming that (L 2 ) w dLL 2 Q(L), we expect the 
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FIG. 5: Boundary effects near the wall determine wall position. In (a), s* eq is determined empirically from MC 
simulations captured in the wall frame and numerically using the finite segment method. In (b), MC simulations are 
plotted in the lab frame. In (c), the occupancy of the last site in the wall frame sn-i = 1 — w+/w- is plotted as a 
function of the mean wall position (N). The parameters a = 1, fc_ = 0.01, k+ = 0, and w+ = 0.001 were used. 



variance of the wall position to be approximately 

y 2 = (t 2 \ - (t) 2 = — 

[ > [ > ew^F'{st q )s' B {{L)-e{m + l)) 

(21) 

In Figure [51 we plot an example distribution Q(L) 
found using both Monte Carlo simulations and from 
(fT8]) . We have aligned the distributions such that 
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FIG. 6: Probability density of the wall position, Q(L), 
is plotted as a function of the deviation from the mean 
wall position (L). The parameters a = 1, w+ = 0.005, 
w- = 0.01, k- = 0.01, and k+ = yield (L) » 350. The 
distribution predicted from (|18|l is a close approximation 
to that derived from MC simulations. 

their maxima coincide. Using (|2ip . the standard 
deviation E rs 0.1215, in good agreement with 
the standard deviation found from MC simulations 
E w 0.143. Since E/(L) - 0.14/3.5 < 1, the wall 
is fairly stable and not likely to fall off the injection 
end of the lattice except on exponentially long time 
scales. 



B. Time-dependent mean wall positions 

In the previous section, we explored the depen- 
dence of the statistically stationary mean wall po- 



sition on the model parameters. However, a stable 
mean wall position may not always exist. In this 
section, we use a FSMFT to determine the stability 
of the wall and the conditions under which a perma- 
nent net wall drift might arise. 

The motion of the wall can be understood com- 
pletely in terms of the outer solution sb(x) - given 
by inverting (fT4")) - and the particle density inside 
finite segment. First, we consider some impor- 
tant properties of Sb(x). Equation (TT¥|) admits two 
branches to the bulk solution, sb(x), because the 
argument of ln[| |] can either be positive or nega- 
tive. The argument approaches zero as s approaches 
St = k + /{k + + k—), the density arising from Lang- 
muir kinetics alone. When we invert x{sb) to find 
sb{x), we see that for increasing x, one branch of 
the density profile Sb{x) approaches sr asymptoti- 
cally from below, and a second branch approaches sr 
asymptotically from above. A representative Sb{x) 
is plotted in Figure [7] Notice that s'(x) > in the 
lower branch, and s'(x) < in the upper branch 
[32| . If a > sr, the steady state density profile will 
lie on the upper branch and the bulk density will 
have values satisfying sr < sb(x) < a. If the in- 
jection rate a < sr, the steady state density profile 
will lie on the lower branch and the bulk density will 
attain values a < sb(x) < sr- 

We are now ready to derive conditions for the 
existence of a fixed mean wall position and stabil- 
ity criteria. In the adiabatic approximation (|20p . 
m is the number of sites in the finite segment and 
e = 1/Nq. This equation expresses the velocity of 
the wall at position eN = L in terms of a particle 
density at position (N — m — l)e. Since s* eq is the 
value of s* that puts no net drift on the wall, i.e., 
w + — W-(l — F(s* q )) = 0, we can expand V(L) from 
([20)1 in a Taylor Series about s* q to find 

V{L) » ew^F'( s ; q )(s*(L - e(m + 1)) - s* eq ). (22) 
Because F is a monotonically increasing function 
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FIG. 7: Four possibilities for the wall motion. Arrows 
indicate the mean wall motion if it is at position i. (a) 
the wall could have no net drift and a stable fixed mean 
position, with small perturbations to the wall position 
decaying over time, (b) the wall could have no net drift, 
but an unstable fixed mean position, with small pertur- 
bations to its position causing it to drift indefinitely to 
the left or right, (c) the wall could drift indefinitely to 
the left and (d) the wall could drift indefinitely to the 
right. The outer solution sb(x) is shown by the dotted 
curves. 



and W- > 0, the wall drifts to the right if s*(L — 
e(m+ 1)) > s* 9 , to the left if s*(£-e(m + l)) < s* eq 
and has a fixed mean position if s*{L — s(m + 1)) = 
s* q . If we now assume that m is sufficiently large 
so that the point L — e(m + 1) lies outside of the 
boundary layer, then s*(L — e(m + 1)) can be well 
approximated by the outer solution given by sb (x) , 
i.e. sb(L — e(m + 1)) w s*(L — e(m + 1)). Further- 
more, we know that sb{x) satisfies a < sb < sr on 
the lower branch and sr < sb < a on the upper one. 
Therefore, we conclude that if s* eq ^ [a, sr], the wall 
can never have a fixed mean position. In particular, 
for all t 



V{L{t)) >0 
< 



if 
if 



s eq < a,s r , 
s* eq > a, s r , 



(23) 



corresponding to an indefinite rightward and left- 
ward drift, respectively (cf. Fig. [7h and d). If a 
fixed mean position does exist, we can understand 
its stability by considering the sign of dV/dL. If this 
quantity is negative (positive), the position is stable 
(unstable). These possibilities are summarized in 
Figure [7] By differentiating (f2"2")l , we have 



^-^ £w _F'( S * eq ) S ' B (L-e(m+l)). (24) 

Hence, if a mean wall position exists at L, a neces- 
sary and sufficient condition for its stability is 
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(25) 



FIG. 8: In (a), we plot the density profile near the wall 
from FSMFT and from MC simulations. We find that 
sjv-i is less than 0.9 = 1 — w+/w~, the value required 
for the wall to have zero drift. In (b), the position of the 
wall N(t) found from MC simulations is plotted with 
the expected N(t) calculated using FSMFT. Parameter 
values were a — 0.01, w+ — 0.001, W- = 0.01, k+ 
:).()!)() 1 . k- = 0.01. 



In particular, when there is no adsorption (fc + = 
0), the bulk solution Sb{x) decreases monotonically 
from the injection site and any mean wall position 
(L) induced by the kinetics will be deterministically 
stable. 

Figures [8] and [9] compare the result ([2D]) with sim- 
ulation data. Figure [8] shows the results of a MC 
simulation in which the wall particle has a mean left- 
ward drift. In Fig. [8^,, the density profile in the wall 
frame found from MC simulation and that predicted 
using FSMFT are shown. Far from the injection site 
sb asymptotes to sr. Thus, when the wall starts at 
a position Lq ^> 1, we assume s*(L — e(m + 1)) = sr 
and the wall's velocity V is independent of its po- 
sition L. Since sjv-i < 1 — we expect, from 
([20]) . that the net drift on the wall will be nega- 
tive. In Fig. Hb, we compare N(t) — L(t)/e found 
from MC simulations with that calculated assuming 
L(t) = Lo + Vt where V is calculated using (f20|) in 
the large L limit. Similarly, Fig. [9^ shows a den- 
sity profile from MC simulation in the case where 
the wall acquires a mean rightward drift. In Fig. [SJs 
both MC simulations and FSMFT show that in the 
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wall frame, the occupancy of the site adjacent to the 
wall is greater than 1 — and we expect a mean 
rightward drift. In Fig. [9};, we see that this is the 
case, and the predicted time course N(t) = L{t)/e 
is compared with N(t) found using MC simulations. 

In contrast to the case of a static mean wall po- 
sition, when the wall has a position-independent ve- 
locity, V, the diffusion constant D of the wall is given 



by D = e (w + 



■W-F(s r ))/2. The probability 



density Q(L,t) describing wall position then follows 



dQ 
dt 



D 



d 2 Q T BQ 



dL 2 



V 



dV 



the solution of which is 



Q(L,t) 



1 



(L-Vt-Lg)^ 



(26) 



(27) 



We now discuss our results in the context of the 
phase transitions [HI, [HI H^] of the interior den- 
sity. When Langmuir kinetics is coupled to a fixed 
domain TASEP with open boundaries, qualitative 
properties of s(x) can change abruptly when adsorp- 
tion/desorption and injection/ejection rates vary. 
For example, an interior boundary layer separating 
regions of low and high density can suddenly disap- 
pear, replaced with a single region of high density as 
the injection rate a is increased. 

Our moving wall TASEP system coupled with 
Langmuir kinetics does not support the phase struc- 
ture seen in [l]| 0j| HH- Because we limit our- 
selves to k+ < k-, we can see from (TT3"f that when 
sb(x) > 0.5, (corresponding with a high density re- 
gion), s' B (x) > 0. From (f2"S")) , the wall cannot have a 
stable equilibrium position within the high density 
region, and we do not find time-independent den- 
sity profiles with low to high density interior shocks 
(a low-high shock), as is observed in [IH, [lfj. In 
fact, references [29[ show that high- low shocks 
are never stationary in an exclusion process with 
Langmuir kinetics. Therefore, interior shocks are 
never stable in our model system. In our problem, 
the presence of a wall that responds to particle dy- 
namics relaxes any shocks in density that may oth- 
erwise occur in the interior, forcing them to the left 
or right boundaries. 



VI. SUMMARY AND CONCLUSIONS 

Our model of an asymmetric exclusion process 
with Langmuir kinetics and a movable right bound- 
ary, and the corresponding results provide a guide to 
understanding biophysical processes in which many 
processing molecular motors push against a load. 
The detachment and attachment rate of the mo- 
tors, as well as the injection rate at the entry site, 



determine the load the motors can support. If a 
static load particle position is reached, we see that 
the mean wall distance from the injection site sat- 
urates upon increasing injection rate a past about 
0.5. The analyses can be used to predict whether 
biological processes such as ribosome movement and 
filopodia/filament extension continues or reaches a 
static configuration. 

Within our model, we found four parameter 
regimes. In the first regime, (sr < s* eq < ot), the wall 
attains a stable equilibrium position for the wall. In 
the second regime, (a < s* eq < Sr), there is an equi- 
librium, but unstable mean wall position. In the 
third and fourth regimes, (s* ^ [a, sr]), the wall 
will always feel a net drift to the to the right and 
left, respectively. In the latter case, the wall will 
fall off the lattice in a time scaling linearly with the 
starting position. When there is a stable equilibrium 
wall position, we can find the mean wall position (N) 
as a function of the particle injection rate a, the ad- 
sorption and desorption rates k±, and the intrinsic 
hopping rates of the wall w± . Determination of (TV) 
requires accurate evaluation of the particle density 
near the wall. Using a hybrid finite segment/mean 
field approach in the reference frame of the fluctuat- 
ing wall, we accurately determine the particle den- 
sity near the wall, and use this to determine the 
wall's steady-state position. 

When there is no steady state wall position, the 
finite segment mean field approach allows us to es- 
timate the steady state velocity of the wall far from 
the injection site. In our analysis, we assumed that 
the particle density has reached steady state, thus 
ignoring the initial particle density profile and wall 
position. Even in regimes where we expect an equi- 
librium wall position at steady state, if the wall is 
initially near the injection site, and the particle den- 
sity is initially very low, we would expect the wall 
to fall off the lattice before reaching its equilibrium 
position. The times to falling off the lattice may be 
treated with extensions of large deviation theory as 
suggested by Fig. [5] 30] . A number of interesting ex- 
tensions of the free boundary problem arise. For ex- 
ample, we expect for certain parameter regimes that 
slow bottleneck sites [6| can attract the fluctuating 
wall. These features and other novel applications to 
biophysical systems deserve investigation. 
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FIG. 9: When sr > sl q and a > sl q , the wall escapes. In (a), far from the injection site, the particle density approaches 
sr = 0.029 as predicted by analytic theory. In (b), we use FSMFT and MC to find sjv-i > 0.9 = 1 — w+/w-, the 
value for which the wall's drift would be zero. In (c), we show N(t) to compare the escape velocity calculated from 
finite segment analysis to the escape velocity found in simulations. Although the value of sjv-i found by FSMFT 
differs from that found in MC simulations by only 0.4%, the calulated velocities differ by 17%. Parameter values 
were a = 0.3, w+ = 0.001, W- = 0.01, k+ = 0.0003, and k- = 0.01. 



VII. APPENDIX 



We can take the continuum limit of (J5]) by denning 
x = (i — l)e where e is the lattice spacing. We find, 

= £«'(2s - 1) + ^" - + Jt + (1 - s) = 0. 

(28) 

I 

The left hand boundary condition, (J6j> becomes 

ds(0,t) 



Of 



a(l - s(0)) - fc-s(0) + fc+(l - s(0)) - s(0)(l - s(e)) = 0. (29) 



When significant changes in the solution near x = vary over a length scale that is > 0(e), this equation is 
well approximated by 

= a(l - s(0)) - ks(0) + fe+(l - «(0)) - s(0)(l - s(0) - es'(0)) = 0. (30) 

The right hand boundary condition analogous to (|10[) . becomes 
ds(L,t) 



dt 



-ks(L) + k+(l - s(L)) + s(L - e)(l - s(L)) + ws{L -s)(l- s(L)) - w+s(L) = 0. (31) 



Again, when significant changes in the solution near x = L vary over a length scale that is > 0(e), this 
equation is well approximated by 

= -k_s(L) + fc+(l - s(L)) + {s(L) - es'(L))(l - s(L))+w_ (s(L) - es'(L))){l - s(L)) - w+s(L) = 0, 

at 

(32) 

at the free boundary L(t), where s(L) defines the particle density at the position just to the left of the wall. 
In the wall frame, the continuum limit of equation ([H]) is 

= eK - (1 + te-(l - sjv-i))]s'(1 - 2s) + j [1 + w+ + w-(l - s N -i)]s" - fc_s + *+(l - s). 

, (331 



In the wall frame, an interior particle shifts to it does with a (normalized) rate of unity, or when 
the right when it either hops to the right, which 
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the wall hops to the left, which it does with rate 
it;_(l — Sjv-i) = w + in steady state. Similarly, a 
particle shifts to the left when it hops to the left 
or when the wall hops to the right, which it does 
with rate w+. In the bulk, where s'(x) = O(l), the 
diffusive term is small and can be neglected. The 
only term we retain that depends on hopping rates 
is — 1 — u>+)(l — 2s)s', which is equal to the 
value of the corresponding term in the lab frame, 
— (1 — 2s)s'. The bulk density is described in both 
frames by 

£s'(l-2a) + ks-k+(l-s)=0. (34) 

In our problem, the second order term in (|28p be- 
comes important in the right hand boundary layer, 
and in the left hand boundary layer when there is 
one. On the left, when s ~ 0.5, one cannot assume 
that s varies slowly. Making the change of variables 
x = £X (£ <C 1), equation (|28|) becomes 

2 

f a' {2a - 1) + ^s" - k-a + k+(l - a) = 0; (35) 

furthermore, we know that e is necessarily on the 
order (fc_ + fc + ). Since the first order term becomes 
very small as s — > 0.5, the second order term must 



match either the adsorption or desorption term. In 
this case, the second order term will be balanced 
when £ ~ yfe. Therefore, we expect a boundary 
layer of width 0{yfe) to arise near the injection site 
if a > 0.5, st < 0.5 or if a < 0.5 and s r > 0.5. 

While the boundary layer on the left can be cap- 
tured using a second order continuum equation, the 
boundary layer on the right cannot. In equation 
(|28j) , we kept terms only up to order e 2 in our expan- 
sion s(x+e) — s(x)+es'(x) + ^-s"(x) + ^rS^ 3S) (x) + .... 
The boundary layer on the right hand side arises 
to join the outer solution with the boundary con- 
dition s(L) = 1 — W+/W-. Making the substitution 
X = x/£, and matching first and second order terms 
in (J3SJ), we find that £ = 0(e). In the boundary layer 
on the left, s « 0.5, and we assume that the term 
e(2s — l)s' is relatively small. When we match the 
second order term, ea"/2, with the adsorption and 
desorption terms, fc+(l — s) and fc_s, we find that 
£ ~ \fe. However, on the right, we cannot assume 
that s w 0.5. We therefore assume that the leading 
terms are e(2s — l)s' and ea"/2, which leads us to 
conclude that the wall-hugging boundary layer has 
width of 0(e). In this case, all terms e n s^ n \X)/n\ 
in the Taylor expansion of (f3"Tj) are 0(1), and con- 
tinuum theory breaks down. 
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Phys. 110, 775 (2003). 
] In [l3l . Il5l ]. an asymmetric exclusion process in a 
fixed domain with open boundaries and Langmuir 
kinetics was studied. Regimes arise in which the po- 
sition of a shock in the particle density becomes in- 
sensitive to the ejection rate /3 on the right, once 
P > 0.5. Because these observations were made 
under the assumption k+ > fc_, and we consider 
k- > k+, by particle hole symmetry, the ejection 
rate in these works corresponds to the in ject ion rate, 
a in our problem. While the authors of [13L Qj| find 



a shock position insensitive to boundary conditions, 
we find an insensitive mean wall position. 
[32] When sb(x) passes though sb(x) = 0.5, the sign of 
s' B (x) changes. The profile s(x) departs from sb(x) 
near sb(x) = 0.5 because one of the assumptions 
used to derive sb{x) - that s' B (x) = 0(1) - becomes 
invalid and the sign of s'(x) fails to change upon 
passing through s(x) — 0.5. See the Appendix for 
further discussion of when the continuum equations 
may cease to be valid. 



